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Gold-palladium (AuPd) nanoparticles have shown significantly enhanced activity relative to monometallic 
Au and Pd catalysts. Knowledge of composition and metal domain distributions is crucial to understanding 
activity and selectivity, but these parameters are difficult to ascertain in catalytic experiments that have 
primarily been devoted to equimolar nanoparticles. Here, we report AuPd nanoparticles of varying Au:Pd 
molar ratios that were prepared by a seed growth method. The selective oxidation of benzyl alcohol was used 
as a model reaction to study catalytic activity and selectivity changes that occurred after varying the 
composition of Pd in bimetallic catalysts. We observed a remarkable increase in catalytic conversion when 
using a 10 : 1 Au:Pd molar ratio. This composition corresponds to the amount of Pd necessary to cover the 
existing Au cores with a monolayer of Pd as a full-shell cluster. The key to increased catalytic activity derives 
from the balance between the number of active sites and the ease of product desorption. According to 
density functional theory calculations, both parameters are extremely sensitive to the Pd content resulting in 
the volcano-like activity observed. 



Optimizing the architectures of bimetallic nanoparticles constitutes a major goal to merge synthesis with 
structure- activity relationships that allow the rational design of catalysts^'^. Gold-palladium nanoparti- 
cles (AuPd NPs) have shown significantly enhanced catalytic activity relative to monometallic Au and 
Pd NP catalysts in prototypical reactions, such as CO oxidation^'^, and in industrially challenging reactions, such 
as vinyl acetate monomer (YAM) synthesis^, the direct synthesis of hydrogen peroxide from H2 and 02^, alkene 
epoxidation^, and the oxidation of hydrocarbons^ and alcohols^'^". Bimetallic nanoparticles are formed by metallic 
domain distributions ranging from homogeneous alloys (AB) to core-shell (Acore@Bsheii or Bcore@Asheii) nano- 
particles. The type of morphology prepared depends on the synthetic protocol, the miscibility of both metals, and 
post- synthesis treatment. Pd and Au are miscible in almost all compositions, forming AuPd alloy NPs^'^^"^^. The 
AuPd alloy NPs most frequently used in catalytic applications contain a 1 : 1 Au:Pd molar ratio. However, it has 
been well documented that the surface of a AuPd alloy NP differs from its corresponding bulk concentration^^'^^. 
Enache et al.^ have shown that Pd segregates towards the surface upon calcination, producing alloy NPs with a Pd- 
rich shell and an Au-rich core. Therefore, the active catalyst has a surface that is significantly enriched with Pd, 
obtained by the metal's migration to the surface during the annealing process. Additionally, core-shell NPs can be 
produced by the reduction of Pd over preformed Au NPs, and vice versa. Henning et aP° have shown the 
preparation of core-shell Au-Pd catalysts by the growth of Pd on icosahedral Au seeds through thermodynamic 
control. In catalytic studies, they identified a 12 atomic-layer-thick Pd shell (2.2 nm) as the optimum catalyst for 
benzyl alcohol oxidation. However, no mechanistic explanation was given for the observed composition effects. 
NPs with different Au:Pd molar ratios have also been suggested as the most active catalysts in other studies^'^^"^\ 
Thus, the question of which is the best Au:Pd ratio in the optimal catalyst for a given reaction and why, still 
remains unanswered. Here we synthesized Aucore@Pdsheii NPs with low Pd loading and investigated how the 
catalytic activity and selectivity were affected both experimentally and theoretically. 

Results and Discussion 

The preparation of bimetallic Au@Pd nanoparticles was based on the reduction of varying amounts of palladium 
over supported gold nanoparticles (Au NPs). The catalyst support chosen for this study is comprised of magnetic 
cores (Fe304 — 10 nm) spherically coated with silica (spheres of —40 nm in diameter)^^. The catalyst support was 
loaded with Au^^ ions (0.9 wt%, determined by ICP OES), and supported Au NPs were formed by controlled 
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Figure 1 | Catalytic performance of the Au@Pd core-shell NP catalysts in 
the oxidation of benzyl alcohol. The amount of Au is fixed (3.4 |amol), 
while the amount of Pd varies from 0 to 40 mol % (i.e., 0 to 1.4 famol). 

reduction in molecular hydrogen. The morphology and size of the 
Au NPs used as seeds for the deposition of Pd were determined using 
a scanning transmission electron microscope (STEM). The Au NPs 
demonstrated a mean diameter of 12.0 ± 3.2 nm, according to rep- 
resentative High Angle Annular Dark Field (HAADF)-STEM images 
and the size distribution histogram fitted to a Gaussian function 
shown in Figure SI. In a second step, Pd (II) acetate solutions in 
benzyl alcohol containing variable amounts of Pd^^ ions were added 
over the Au cores and the reduction of Pd occurred in molecular 
hydrogen for 30 minutes at 100°C and 4 bar H2. The amount of Au 
(in mols) in the resultant Au@Pd NPs remained fixed, while the Pd 
composition varied from 1 to 40 mol % relative to Au, corresponding 
to 0 to 4 monolayers (ML) of Pd. The aerobic oxidation of benzyl 
alcohol was used as a model reaction to study the activity and select- 
ivity of the Au@Pd core-shell NP catalysts (Figure 1). Conversion of 
benzyl alcohol and selectivity obtained by selected Au@Pd core- shell 
NP catalysts and monometallic Au and Pd NP catalysts are presented 
in Table 1. The monometallic Au catalyst promoted high selectivity 
for benzaldehyde formation, but less than 2% conversion (Table 1 
entry 1). Both conversion and selectivity were poor for the mono- 
metallic Pd catalyst prepared by adding the corresponding amount of 
Pd to the support without gold (Table 1 entry 9). The bimetallic 
catalysts exhibited a strong dependence between the activity and 
the amount of palladium added as shell metal. The highest benzyl 
alcohol conversion was achieved by adding 10 mol % Pd in respect to 
the existing amount of Au (i.e., Au:Pd = 10:1) (Table 1, entry 5). A 
sharp increase in benzyl alcohol conversion was observed by adding 



Table 1 | Catalytic oxidation of benzyl alcohol by Au, Pd and Au@ 
Pd catalysts 

Catalyst Selectivity (%) 



Entry 


\xmo\ Au 


[imol Pd 


Conv. (%) 


1 


II 


III 


IV 


1 


3.4 


0 


1.7 




97.6 


0 


2.4 


2 


3.4 


0.039 


14.0 




95.2 


1.7 


3.0 


3 


3.4 


0.20 


41.2 




89.8 


7.6 


2.6 


4 


3.4 


0.26 


67.0 


0.2 


88.2 


7.8 


3.7 


5 


3.4 


0.34 


87.7 


0.4 


86.4 


9.4 


3.9 


6 


3.4 


0.51 


64.6 


0.3 


93.3 


3.0 


3.3 


7 


3.4 


0.78 


39.2 


0.3 


91.8 


4.3 


3.6 


8 


3.4 


1.4 


17.5 


0.6 


95.0 


1.9 


2.5 


9 


0 


3.4 


10.0 


26.7 


67.8 


2.5 


3.1 


10 


0 


6.8 


30.2 


2.3 


89.7 


4.1 


3.9 



l-Toluene; ll-Benzaldehyde; III- Benzoic acid; IV-Benzyl benzoate. Reaction conditions: benzyl 
alcohol (9.6 mmol), T = 1 00°C, p02 = 6 bar, 2.5 h reaction. 



Table 2 | Catalyticoxidationof benzyl alcohol by recycling Au@Pd 
(Au:Pd 10: 1) catalyst 

Selectivity (%) 



Recycle 


Conv. (%) 


1 


II 


III 


IV 


1 


60.1 


0.2 


93.7 


4.5 


1.6 


2 


56.3 


0.3 


96.3 


2.1 


1.3 


3 


65.1 


4.2 


94.3 


0.9 


0.5 


4 


37.7 


2.7 


93.8 


2.4 


1.1 


5 


45.2 


2.2 


93.7 


3.1 


1.0 


6 


50.3 


2.0 


92.5 


4.6 


0.9 



l-Toluene; ll-Benzaldehyde; III- Benzoic acid; IV-Benzyl benzoate. Reaction conditions: benzyl 
alcohol (9.6 mmol), bimetallic catalyst (3.4 \imo\ Au + 0.34 nmolPd), T = 1 00°C, p02 = 6 ba 
2.5 h reaction. 



Up to 10 mol % Pd (or 1 ML Pd), which was followed by a decrease as 
the Pd content was increased, as shown in Figure 1. The selectivity to 
benzaldehyde was only slightly affected by the Pd content of Au@Pd 
catalysts. After each reaction, the catalysts were easily separated by 
placing a magnet on the external reactor wall due to the magnetic 
properties of the catalyst support. The most active catalyst (Au:Pd = 
10:1) was recycled and reused in successive reactions, maintaining 
the high selectivity to benzaldehyde (Table 2). Fluctuation in the 
reaction conversion was reported before for an alloyed AuPd cata- 
lyst^^ and can be related to the presence of amino-groups on the 
catalysts support (used to improve the impregnation of gold on 
the support^^) and possible interactions between those groups and 
the reaction products. 

The metal composition of the most active catalyst was determined 
by Flame Atomic Absorption Spectrometry (FAAS) to consist of 
mole fraction percentages of 89.9% Au and 9.1% Pd (i.e., Au:Pd = 
10:1), which agrees with the nominal composition. This was also the 
case for other catalysts, as shown in Table 3, which suggests depos- 
ition of the desired amount of palladium. 

The most active catalyst was also investigated by X-ray energy 
dispersive spectroscopy in a scanning transmission electron micro- 
scope (XEDS-STEM) to determine both the compositional distri- 
bution between particles and the existence of compositional rich 
domains within each particle. The mean composition (mole fraction) 
of the catalyst particles was determined to be 10.2 ± 6.7% Pd (see ESI 
for the details. Figure S2). The distribution of the metallic compo- 
nents within each particle was investigated using X-ray energy dis- 
persive spectroscopy spectrum imaging (XEDS-SI). The XEDS-SI 
technique enables direct visualization of the nanoscale structural 
and compositional features of the synthesized core-shell particles^^. 
Figure 2 shows a HAADF-STEM image of a quasi- spherical catalyst 
particle with the respective Au and Pd maps. In the compositional 
maps, the bright intensity in each pixel is a measurement of the 
number of X-ray peak counts at each corresponding position. 
Although the number of X-ray counts for Au Loc is larger than for 
Pd La in the investigated material (Figure S2), the brightness distri- 
bution was normalized in the maps shown in Figure 2. In the Au map. 
X-ray counts are measured throughout the particle, resulting in a 
hemisphere- shaped line scan. A different situation is observed in the 
Pd distribution because the number of counts increases steeply at the 
borders of the particle and remains nearly constant throughout 
the remainder of the particle. This distribution shows that the part- 
icles investigated are morphologically structured with an Au-rich 
core and a Pd-rich shelP^. The less active catalyst with Au:Pd = 
10:4 is also morphologically structured with an Au-rich core and 
a Pd-rich shell (Figure S3). Monometallic nanoparticles were not 
found in the investigated samples. 

In the catalytic studies, we observed a volcano -like behavior char- 
acterized by a steady increase in benzyl alcohol conversion as the 
quantity of added Pd (shell metal) was increased relative to Au (core 
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Table 3 Composition of Au and Pd 


in selected Aucore@Pdshel 


catalysts determined by FAAS 






Cotalyst (AuiPcl molor rotio) 


%Au 


%Pd 


Au (jirnol g ) 


ru ^)j,moi g ) 


Au;Pd molor rotio 


Au@Pd(10:l) 


1.42 


0.08 


72 


7.5 


1 :0.10 


Au@Pd(10:1.5) 


1.14 


0.12 


58 


1 1 


1 :0.18 


Au@Pd(10:3) 


0.99 


0.16 


50 


15 


1 :0.30 


Au@Pd (10:4) 


1.28 


0.27 


64 


25 


1 :0.39 



metal), with a consistent decrease in conversion for more than 
10 mol % Pd. If the Pd atoms are primarily on the surface of Au, 
as observed in the morphological studies, the composition of the 
most active catalyst (89.9% Au and 9.1% Pd) correlates with the 
composition required for the complete coverage of Au cores (d = 
12.0 ± 3.2 nm) with one atomic layer of Pd (see Table SI). A special 
behaviour related to the monolayer regime was reported previously 
for other alloys^^"^^. 

Theoretical simulations are fundamental to the understanding of 
the volcano-like behavior reported in Figure 1. Density functional 
theory-based investigations of AuPd systems have concentrated on 
the structure^^"^^, oxygen activation^^, or VAM synthesis^^"^^'^^. Here 
we modelled the reaction process on three representative surfaces 
with an ordered surface alloy (ML denotes the number of Pd mono- 
layers grown on the Au seeds): Au@Pd3Au, Au@Pd(l ML) and Au@ 
Pd(2 ML) (Figure 3); and Pd. 

With these material models and the PBE functionaP^ we inspected 
the paths leading to the formation of the aldehyde, considering the 
following elementary steps: 



Figure 2 | HAADF-STEM image ofa supported catalyst particle Au:Pd = 
10:1 and the corresponding Au and Pd maps. The particle composition 
distribution is observed in the line scans, measured from the regions 
delimited by the lines indicated in both maps. 



(1) Alcohol adsorption: C6H5CH20H(g) + * CgHsCHsOH* 

(2) Oxygen adsorption: 02(g) + 2* ^ 20* 

(3) First H transfer: CgHsCHsOH* + O* ^ CgHsCHsO* + OH* 

(4) Second H transfer: C6H5CH2O* + OH* CgHsCHO* + 
H2O* 

(5) Water desorption: H2O* H20(g) + * 

(6) Aldehyde desorption: CgHsCHO* C6H5CHO(g) + * 

(7) Unselective path: CgHsCHO* + H2O* CgHsCHOH* + 
OH* 

The reaction profile is presented in Figure 4 for the bimetallic 
systems only. The activation barriers, not shown in the figure, and 
all data for the Pd-only system, are presented in Table S2. Analysis of 
the results in Figure 4 leads to the conclusion that adsorption of the 
alcohol is difficult (endothermic by 0.45 eV) on Au@Pd3Au, and 
thus the reaction cannot progress adequately in such structures 
because the number of molecules adsorbed on the surface is small. 
Increasing the amount of Pd on the surface improves adsorption 
rendering it exothermic ( — 0.31 eV) as observed for Au@ 
Pd(l ML). When additional Pd is present in the slab, the binding 
energy is very strong, — 1.69 eV for Au@Pd(2 ML). The reaction has 
easily attainable energy barriers for all of the steps (see Table S2). 
Water desorption (step 5) is exothermic for systems with Pd shells up 
to one monolayer ( — 0.66 and —0.22 eV for Au@Pd3Au and Au@ 
Pd(l ML), respectively), and endothermic for NPs having higher Pd 
content. In the reaction mechanism, the OH intermediate is a short 
living species due to the difficulties found in its formation and the 
easiness at which the next reaction takes place. Therefore, the main 
path for the unselective process is the back transfer of one proton 
from water, i.e. step 7. The reduced energy barrier for water desorp- 
tion guarantees the selectivity towards benzaldehyde because the 
barrier for this step is, in all cases, lower than H transfer from water 
to the aldehyde (Table S2). The final aldehyde product is easily 
desorbed from the Au@Pd3Au model, but its binding energy is large 
on Pd monolayers: 0.8 and 1.1 eV, for 1 and 2 ML, respectively. As 
for the activity, reactant adsorption (1.7 vs. 0.3 eV) and product 
desorption (1.1 vs. 0.8 eV) are far two strong on the two monolayer 
system compared to the single layer to ensure a rapid recycling of the 
clean surface. 

The computational approach is in line with the mechanistic 
investigations in the literature. For instance, Enache et al.^ found zero 
order for oxygen, which is in agreement with the ease of dissociation 
of oxygen reported in previous calculations^^. In addition, the appar- 
ent activation energy reported was 0.47 eV. This value is similar to 
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Figure 3 | The three model systems employed in the DFT calculations. 

Left a PdsAu overlayer on Au; center a Pd monolayer on Au and right two 
Pd monolayers on Au. Blue spheres stand for Pd atoms and yellow (olive) 
for Au. 
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Figure 4 | Reaction profile for the selective oxidation of benzyl alcohol to 
the corresponding aldehyde. The insets schematically represent the 
different steps in the reaction for the Au@Pd( 1 ML) model: blue spheres 
Pd, yellow Au, grey C, red O, white H (see also Figure S4). 

the AG for the desorption of benzylaldehyde in Au@Pd(l ML) once 
dispersion and entropic terms are considered (0.46 eV). Pd grown on 
Au has an upward shift of the Pd band by +0.30 eV^^. Strain induced 
by Au on Pd should increase the reactivity as the Au lattice is larger 
than that of Pd. Indeed enhanced alcohol adsorption is found for the 
two Au@Pd(2 ML) model, -L69 vs. -0.86 eV, see Table S2. 
However, the strong Pd-Au ligand effect is observed and indeed 
the ligand contribution affects alcohol adsorption to such an exten- 
sion that the binding is —0.31 eV only. The calculations including 
entropy and dispersion were done as follows: AG = AE — TAS; 
where S was only considered for the gas-phase species and AE 
included the dispersion terms calculated through the dispersion 
model introduced by Grimme et al.^^ and the metal parameters from 
Tkatchenko et aP^. In the case of Pd, to analyze unselectivity, water 
desorption and unselective H transfer need to be considered. If entro- 
pic and dispersion contributions are considered desorption is endo- 
thermic by 0.3 eV while the barrier for proton transfer is as low as 
0.14 eV, therefore selectivity is compromised. 

The analysis of the reaction network accounts for the observation 
of a maximum in the experimental catalytic results. On the left side of 
the volcano-like plot (Figure 1) corresponding to low Pd loadings, 
DFT results indicate that adsorption of the alcohol is controlling the 
activity. The electronic properties of AuPd alloys have been exten- 
sively analyzed in the literature^^"^^. Actually, only sites with no Au 
would be active for benzyl alcohol adsorption, and thus, activity is 
proportional to the number of Pd ensembles (i.e., at least 6 Pd atoms 
on the surface) that can accommodate benzyl alcohol. However, 
increasing benzyl alcohol adsorption subsequently increases the 
desorption energy required to clean the surface from the aldehyde 
product. This is a drawback when having more than 1 ML of Pd, as 
the binding energy of the aldehyde is far too large to remove the 
product from the surface and leave clean ensembles for the next 
catalytic cycle. As a consequence, the best balance is obtained when 
1 ML of Pd covers the Au nanoparticles because benzyl alcohol 
adsorption can then occur everywhere and desorption is still feasible. 
For the Pd-only system, the barrier for the second H transfer is the 
most energy- demanding step, 1.32 eV (Table S2), and thus implies 
that the reaction is slow at the temperatures at which AuPd alloys are 
active. In addition, once the aldehyde is formed, the water desorption 
energy and the barrier for H stripping are comparable (0.10 vs. 
0.14 eV, respectively), and therefore, compromising selectivity 



determining step is not as favorable as for the alloys. When extending 
the study to other unsaturated alcohols the selectivity controlled by 
the competition between steps (5) and (7) would follow the same 
general rule. As for the activity, the adsorption of these unsaturated 
alcohols (like crotyl alcohol) needs smaller ensembles than those of 
the corresponding benzyl alcohol, but again the number of sites 
would maximize for a composition of 1 ML of Pd on top of the Au 
core. Therefore, the present results can be extrapolated to other 
unsaturated alcohols. 

Conclusions 

We have shown that the rational design based on the density func- 
tional theory and the controlled synthesis of Au@Pd core-shell cat- 
alysts with specific compositions and metal domain distributions can 
be used as a strategy to enhance the catalytic properties of Au while 
maintaining its high selectivity for the oxidation of benzyl alcohol to 
benzaldehyde. This finding shows that many opportunities remain to 
obtain enhanced properties in bimetallic nanomaterials by tailoring 
the amount of each metal and the metal domain distribution, and 
opens new possibilities in the field of nanomaterials synthesis and 
design. 

Methods 

Synthesis of the catalyst support. The core-shell Fe304@Si02 nanocomposite was 
prepared and modified with (3-aminopropyl)triethoxysilane as reported in the 
literature^*^. 

Synthesis of the supported Au monometallic catalyst. The catalyst support, silica- 
coated magnetite nanoparticles previously modified with (3- 

aminopropyl)triethoxysilane (500 mg), was added to an aqueous solution of HAuCU 
(40 mL, 0.215 g L"^) with pH 6. The mixture was stirred at 25°C for 2 h. The solid 
was then collected magnetically from the solution and washed twice with hot distilled 
water and twice with ethanol. The reduction step was performed in a Fischer-Porter 
glass reactor pressurized with hydrogen gas. The gold precursor was dispersed in 
cyclohexane and heated to 100°C under H2 at 4 bar. The mixture was stirred for 2 h, 
and the color of the solid changed from brown to dark red. 

Synthesis of the supported AuPd bimetallic catalysts. The preparation of bimetallic 
catalysts consisted of an addition of the desired amount of a Pd(OAc)2 solution in 
benzyl alcohol (20 to 800 [iL, [Pd'+] = 1.68 X 10"' mol L'^) to 75 mg of Au 
monometallic catalyst (3.4 |amol Au), which was then added to a Fischer-Porter 
100 mL glass reactor. Benzyl alcohol was added to achieve a total volume of 1.0 mL 
and the reactor was pressurized with hydrogen gas at 4 bar. The mixture was stirred 
for 30 min at 100°C. The color of the solid changed to black. The reactor was cooled to 
room temperature, the hydrogen gas was relieved, and the reaction mixture was kept 
inside the reactor for catalytic studies. 

Catalytic experiments. All reactions were performed using a modified Fischer- 
Porter 100 mL glass reactor. In a typical reaction, the glass reactor was loaded with the 
supported Au catalyst (75 mg, 3.4 |amol Au) and benzyl alcohol (9.6 mmol), or 
the bimetallic Au-Pd catalyst-containing mixture was prepared as described above. 
The reactor was purged three times with O2, leaving the vessel at 6 bar. The 
temperature was maintained at 100°C with an oil bath on a hot stirring plate 
connected to a digital controller (ETS-D5 IKA). The reactions were constantly stirred 
using Teflon-coated magnetic stir bars for the desired amount of time. The catalyst 
was recovered magnetically by affixing a magnet onto the reactor wall. The products 
were collected with a syringe and analyzed by gas chromatography (GC) and gas 
chromatography- mass spectrometry (GC-MS). 

Morphological investigation. Catalyst samples were prepared by sonicating the 
catalyst powder in isopropanol. A drop of the resulting dispersion was placed on a 
thin carbon film, which was deposited on standard 400-mesh TEM copper grids and 
air-dried. The size distribution of the supported Au core particles was investigated by 
HAADF-STEM using a JEOL-JEM 2100F TEM microscope available at the LNNano 
Laboratory (CNPEM, Campinas, Brazil). The mean diameter distribution was 
measured by HAADF image analysis using Image-Pro Plus 6.0 software. Two 
hundred nanoparticles were measured to build the size distribution histogram that is 
presented. Compositional distribution was measured on the same instrument by 
acquiring individual XEDS-STEM spectra of a sample of 62 supported particles. 
Particle compositional maps were acquired by XED-Spectrum Imaging with a Digital 
Micrograph 1.8 system (Gatan Inc.) controlling a Thermo-Noran XEDS. The Sis were 
acquired at a 15 X 15 pixel resolution with a dwell time of 2 s/pixel, using the drift 
correction facility every 100 s. XED spectra were acquired in the 0-10 keV energy 
range at each pixel. Electron probe sizes of approximately 0.7 nm were obtained by 
operating the JEM 21 OOF instrument in STEM mode, allowing enough current 
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density at each sample point to acquire statistically significant X-ray counts for the 
investigated elements. 

Computational Methods. Density Functional Theory was applied to slabs modeling 
the system. The VASP program was employed to this end^^'^°. The functional of 
choice was PBE^^. Inner electrons were replaced by frozen cores PAW*^ and valence 
electrons were expanded in plane waves with cut-off energy of 450 eV. As the 
particles are large enough only (111) surfaces, the most stable for fee metals, have been 
employed in the simulations. Three different models were employed to represent the 
surface of the Pd grown on Au particles: the first one is a partial monolayer on the 
surface of Au@Pd3Au composition, this is following one of the ordered alloys, a single 
Pd monolayer on Au, and a Pd bilayer on Au (see Fig. 3) and Pd alone. The slabs 
contain four layers of metal in total and in all cases overlayers follow the lattice 
parameter of Au. The vacuum space interleaving the slabs is at 15 A. As the adsorbates 
are large we have employed a surface reconstructions p(6 X 4) that can account for the 
adsorption of the benzyl alcohol and its derivates. On these supercells the k-point 
sampling employed are 2 X 5 X 1 Monkhorst-Pack points centered at gamma point*^. 
Transition state searches were performed with the climbing image version of the 
Nudged Elastic Band method*^. 
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